arXiv:hep-ex/0503030vl 18 Mar 2005 


Observation of ppn 0 and ppi 7 in i/j' decays 

M. Ablikim 1 , J. Z. Bai 1 , Y. Ban 11 , J. G. Bian 1 , X. Cai 1 , J. F. Chang 1 , H. F. Chen 17 , 

H. S. Chen 1 , H. X. Chen 1 , J. C. Chen 1 , Jin Chen 1 , Jun Chen 7 , M. L. Chen 1 , 

Y. B. Chen 1 , S. P. Chi 2 , Y. P. Chu 1 , X. Z. Cui 1 , H. L. Dai 1 , Y. S. Dai 19 , Z. Y. Deng 1 , 

L. Y. Dong la , Q. F. Dong 15 , S. X. Du 1 , Z. Z. Du 1 , J. Fang 1 , S. S. Fang 2 , C. D. Fu 1 , 

H. Y. Fu 1 , C. S. Gao 1 , Y. N. Gao 15 , M. Y. Gong 1 , W. X. Gong 1 , S. D. Gu 1 , Y. N. Guo 1 , 
Y. Q. Guo 1 , Z. J. Guo 16 , F. A. Harris 16 , K. L. He 1 , M. He 12 , X. He 1 , Y. K. Heng 1 , 

H. M. Hu 1 , T. Hu 1 , G. S. Huang 16 , X. P. Huang 1 , X. T. Huang 12 , X. B. Ji 1 , C. H. Jiang 1 , 

X. S. Jiang 1 , D. P. Jin 1 , S. Jin 1 , Y. Jin 1 , Yi Jin 1 , Y. F. Lai 1 , F. Li 1 , G. Li 2 , H. H. Li 1 , 

J. Li 1 , J. C. Li 1 , Q. J. Li 1 , R. Y. Li 1 , S. M. Li 1 , W. D. Li 1 , W. G. Li 1 , X. L. Li 8 , X. Q. Li 10 , 

Y. L. Li 4 , Y. F. Liang 14 , H. B. Liao 6 , C. X. Liu 1 , F. Liu 6 , Fang Liu 17 , H. H. Liu 1 , 

H. M. Liu 1 , J. Liu 11 , J. B. Liu 1 , J. P. Liu 18 , R. G. Liu 1 , Z. A. Liu 1 , Z. X. Liu 1 , F. Lu 1 , 

G. R, Lu 5 , H. J. Lu 17 , J. G. Lu 1 , C. L. Luo 9 , L. X. Luo 4 , X. L. Luo 1 , F. C. Ma 8 , 

H. L. Ma 1 , J. M. Ma 1 , L. L. Ma 1 , Q. M. Ma 1 , X. B. Ma 5 , X. Y. Ma 1 , Z. P. Mao 1 , 

X. H. Mo 1 , J. Nie 1 , Z. D. Nie 1 , S. L. Olsen 16 , H. P. Peng 17 , N. D. Qi 1 , C. D. Qian 13 , 

H. Qin 9 , J. F. Qiu 1 , Z. Y. Ren 1 , G. Rong 1 , L. Y. Shan 1 , L. Shang 1 , D. L. Shen 1 , 

X. Y. Shen 1 , H. Y. Sheng 1 , F. Shi 1 , X. Shi llc , H. S. Sun 1 , J. F. Sun 1 , S. S. Sun 1 , 

Y. Z. Sun 1 , Z. J. Sun 1 , X. Tang 1 , N. Tao 17 , Y. R„ Tian 15 , G. L. Tong 1 , G. S. Varner 16 , 
D. Y. Wang 1 , J. Z. Wang 1 , K. Wang 17 , L. Wang 1 , L. S. Wang 1 , M. Wang 1 , P. Wang 1 , 

P. L. Wang 1 , S. Z. Wang 1 , W. F. Wang 1 '', Y. F. Wang 1 , Z. Wang 1 , Z. Y. Wang 1 , 

Zhe Wang 1 , Zheng Wang 2 , C. L. Wei 1 , D. H. Wei 1 , N. Wu\ Y. M. Wu 1 , X. M. Xia 1 , 

X. X. Xie 1 , B. Xin 86 , G. F. Xu 1 , H. Xu 1 , S. T. Xue 1 , M. L. Yan 17 , F. Yang 10 , H. X. Yang 1 , 
J. Yang 17 , Y. X. Yang 3 , M. Ye 1 , M. H. Ye 2 , Y. X. Ye 17 , L. H. Yi 7 , Z. Y. Yi 1 , C. S. Yu 1 , 
G. W. Yu 1 , C. Z. Yuan 1 , J. M. Yuan 1 , Y. Yuan 1 , S. L. Zang 1 , Y. Zeng 7 , Yu Zeng 1 , 

B. X. Zhang 1 , B. Y. Zhang 1 , C. C. Zhang 1 , D. H. Zhang 1 , H. Y. Zhang 1 , J. Zhang 1 , 

J. W. Zhang 1 , J. Y. Zhang 1 , Q. J. Zhang 1 , S. Q. Zhang 1 , X. M. Zhang 1 , X. Y. Zhang 12 , 

Y. Y. Zhang 1 , Yiyun Zhang 14 , Z. P. Zhang 17 , Z. Q. Zhang 5 , D. X. Zhao 1 , J. B. Zhao 1 , 

J. W. Zhao 1 , M. G. Zhao 10 , P. P. Zhao 1 , W. R, Zhao 1 , X. J. Zhao 1 , Y. B. Zhao 1 , 

Z. G. Zhao le , H. Q. Zheng 11 , J. P. Zheng 1 , L. S. Zheng 1 , Z. P. Zheng 1 , X. C. Zhong 1 , 

B. Q. Zhou 1 , G. M. Zhou 1 , L. Zhou 1 , N. F. Zhou 1 , K. J. Zhu 1 , Q. M. Zhu 1 , Y. C. Zhu 1 , 

Y. S. Zhu 1 , Yingchun ZlnRQ Z. A. Zhu 1 , B. A. Zhuang 1 , X. A. Zliuang 1 , B. S. Zou 1 

(BES Collaboration) 

1 Institute of High Energy Physics, Beijing 100049, People’s Republic of China 
2 China Center for Advanced Science and Technology 
(CCAST), Beijing 100080, People’s Republic of China 
3 Guangxi Normal University, Guilin 541004, People’s Republic of China 
4 Guangxi University, Nanning 530004, People’s Republic of China 
5 Henan Normal University, Xinxiang 453002, People’s Republic of China 
6 Huazhong Normal University, Wuhan 430079, People’s Republic of China 
' Hunan University, Changsha 410082, People’s Republic of China 
8 Liaoning University, Shenyang 110036, People’s Republic of China 
9 Nanjing Normal University, Nanjing 210097, People’s Republic of China 
10 Nankai University, Tianjin 300071, People’s Republic of China 


1 


11 Peking University, Beijing 100871, People’s Republic of China 

12 Shandong University, Jinan 250100, People’s Republic of China 

13 Shanghai Jiaotong University, Shanghai 200030, People’s Republic of China 

14 Sichuan University, Chengdu 610064, People’s Republic of China 

15 Tsinghua University, Beijing 100084, People’s Republic of China 

16 University of Hawaii, Honolulu, HI 96822, USA 
17 University of Science and Technology of China, Hefei 230026, People’s Republic of China 
18 Wuhan University, Wuhan 430072, People’s Republic of China 
19 Zhejiang University, Hangzhou 310028, People’s Republic of China 
a Current address: Iowa State University, Ames, IA 50011-3160, USA. 
b Current address: Purdue University, West Lafayette, IN 47907, USA. 
c Current address: Cornell University, Ithaca, NY 14853, USA. 
d Current address: Laboratoire de l’Acceleratear Lineaire, F-91898 Orsay, France. 
e Current address: University of Michigan, Ann Arbor, MI 48109, USA. 

* Current address: DESY, D-22607, Hamburg, Germany. 

(Dated: February 7, 2008) 

Abstract 

The processes A —» PP' 1 r° and ip' —> ppr) are studied using a sample of 14 x 10 6 ip 1 decays collected 
with the Beijing Spectrometer at the Beijing Electron-Positron Collider. The branching fraction 
of ip' —► ppTT 0 is measured with improved precision as (13.2 ± 1.0 ± 1.5) x 10~ 5 , and ip' —> ppr) is 
observed for the first time with a branching fraction of (5.8 ± 1.1 ± 0.7) x 10 -5 , where the first 
errors are statistical and the second ones are systematic. 

PACS numbers: 13.25.Gv, 12.38.Qk, 14.20.Gk, 14.40.Cs 
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I. INTRODUCTION 


There are some long-standing puzzles in the decays of vector charmonia, in particular the 
“p7r puzzle” between ip' and J/ip decays and the possible large charmless decay branching 
fraction of the ip". Following the suggestion in Ref. [1] that the small ip' —> pn branching 
fraction is due to the cancellation of the S- and D- wave matrix elements in ip' decays, it was 
pointed out that all ip' decay channels should be affected by the same S- and D- wave mixing 
scheme, and thus, in general, the ratios between the branching fractions of ip' and J/ip decay 
into the same final states may have values different from the “12% rule”, expected for pure 
IS and 2 S states [2], This scenario also predicts ip" decay branching fractions since ip" 
would also be a mixture of S- and D- wave charmonia, and it was suggested that many J/ip 
and ip', as well as ip", decays should be measured to test this. For the channels that have 
been measured, ip' decays are found to be either suppressed (like vector-pseudoscalar, vector- 
tensor), or enhanced (like or obey the 12% rule (like baryon-antibaryon pairs). In 

this paper, we analyze three-body decays of ip' into a pp pair and a 7r° or rp 

The J/ip and ip' decays into ppn° and pprj are expected to be dominated by two-body 
decays involving excited nucleon states. These states play an important role in the un¬ 
derstanding of nonperturbative QCD [3-6]. However, our knowledge on N* resonances, 
based on nN and jN experiments [7], is still very limited and imprecise. Studies of N* reso¬ 
nances have also been performed using J/ip events collected at the Beijing Electron-Positron 
Collider (BEPC) [8-10], which provides a new method for probing N* physics. Recently, 
based on 58 million J/ip events collected by BEijing Spectrometer (BESII), a new N* peak 
with a mass at around 2065 MeV/c 2 was observed [10]. This may be one of the “missing 
N* states” around 2 GeV/c 2 that have been predicted by the quark model in many of its 
forms [4, 11, 12]. However, due to its large mass, the production of this IV* (2065) in J/ip 
decays is rather limited in phase space, and a similar search for it in ip' decays, which has a 
larger phase space available may be helpful, although the production rate may be small due 
to the large decay rates for ip' into final states with charmonium. 

In a recent paper [13], it is predicted that pp in ip' — > pp-ix 0 may form iso-vector bound 
states near threshold. These states can also be searched for. 

Experimentally, ip' —> ppir° was studied by Mark-II in 1983, with only 9 events ob¬ 
served [14], and the branching ratio was found to be (1.4 ± 0.5) x 10 -4 . ip' — > ppp has not 
been observed before. 


II. DETECTOR AND DATA SAMPLES 


The data used in this analysis were taken with the BESII detector at the BEPC storage 
ring at a center-of-mass energy corresponding to AT,//. The data sample corresponds to a 
total of (14.0 ± 0.6) x 10 6 ip' decays, as determined from inclusive hadronic events [15]. 

BES is a conventional solcnoidal magnet detector that is described in detail in 
Refs. [16, 17]. A 12-layer vertex chamber (VTC) surrounding the beam pipe provides trigger 
information. A forty-layer main drift chamber (MDC), located radially outside the VTC, 
provides trajectory and energy loss ( dE/dx ) information for charged tracks over 85% of the 
total solid angle. The momentum resolution is cr p /p = 0.017a/1 + p 2 (p in GeV/c), and the 
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dE/dx resolution for hadron tracks is ~ 8%. An array of 48 scintillation counters surround¬ 
ing the MDC measures the time-of-flight (TOF) of charged tracks with a resolution of ~ 200 
ps for hadrons. Radially outside the TOF system is a 12 radiation length, lead-gas barrel 
shower counter (BSC). This measures the energies of electrons and photons over ~ 80% of 
the total solid angle with an energy resolution of oe/E = 22 %/y/E (E in GeV). Outside 
of the solenoidal coil, which provides a 0.4 Tesla magnetic field over the tracking volume, is 
an iron flux return that is instrumented with three double layers of counters that identify 
muons of momentum greater than 0.5 GeV jc. 

Monte Carlo (MC) simulation is used for the determination of the invariant mass res¬ 
olution and detection efficiency, as well as the study of background. The simulation of 
the BESI1 detector is Geant3 based, where the interactions of particles with the detector 
material are simulated. Reasonable agreement between data and Monte Carlo simulation 
is observed [18] in various channels tested including e + e~ —> (7)e + e _ , e + e~ —> (7 )p + p~, 
J/i/> —> pp and 4>' —> 7r + 7r~ J/x/>, J/ij; —■> £ + £~ (£ = e, /1). 

The signal channels ip' —> ppn°, n° —> 27 and ip' —» ppp, p —> 27 are generated with 
a phase space generator, giving similar pn, pn, pp, and pp mass distributions to those 
observed in data. For 7r° and p channels, 100 000 events each are simulated. To study 
possible background in our analysis, 20 000 ip' —> n°n°J/ip, J/ip —> pp, 20 000 ip' —> 7T 0 n°pp, 
and 30 000 ip' —> 7y c j, Xcj —■► 7 J/^, J /ip —»• pp (J = 0,1, 2) events are generated. 

III. EVENT SELECTION 

The final states in which we are interested contain two photons and two charged tracks. 
The number of charged tracks is required to be two with net charge zero. Each track should 
have good quality in track fitting and satisfy | cos#! < 0.8, where 0 is the polar angle of the 
track measured by the MDC. 

A neutral cluster in the BSC is considered to be a photon candidate when the angle 
between the nearest charged track and the cluster in the xy plane is greater than 15°, the 
first hit appears in the first five layers of the BSC (about six radiation lengths of material), 
and the angle between the cluster development direction in the BSC and the photon emission 
direction in xy plane is less than 37°. Two or three photon candidates are allowed in an 
event, but the two with the largest energies are selected as 7r° or p decay candidates, and 
both of them are required to have energy greater than 60 MeV. 

A likelihood method is used for discriminating pion, kaon, proton, and antiproton tracks. 
For every charged track, we define an estimator as E l = P l / ]>V P l , where P l is the prob¬ 
ability under the hypothesis of being type i, i = 3,4,5 for 7r, K and p or p hypotheses, 
respectively, and P l = Y[jPj( x j)- Here Pj is the probability density for the hypothesis of 
type i, associated to the discriminating variable Xj. Discriminating variables used for each 
charged track are time of flight in the TOF (TOF-T) and the pulse height in the MDC 
(dE/dx). By definition, pion, kaon, proton, and antiproton tracks have corresponding E l 
values near one. In this analysis, both tracks are required to have E 5 > 0.6. 

A four-constraint (four momentum conservation) kinematic fit is made with the two 
charged tracks and two photon candidates; the confidence level of the y 2 fit is required to be 
greater than 1%. A similar fit assuming the two charged tracks are K + K~ is also performed, 
and the y 2 of ■?// —> 77 pp should be smaller than that of i/)' —» 77 K + K~. 
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The scatter plot of the pp invariant mass versus that of the two photon candidates for 
events satisfying the above selection criteria is shown in Fig. 1(a). The two bands with 
m 77 values near 0.135 GeV/c 2 and 0.547 GeV/c 2 are ip' —> pp7r°, 7r° —> 2y and ip' —> PWh 
p —> 2 7 candidates, respectively. The band corresponding to pp mass around 3.1 GeV/c 2 
is from ip' -> 7y cJ , y cJ -> yJ/^, J/ip ^ pp (J = 0,1,2), and ip' -> //J/-0, J/ip -> 
pp. The broadly distributed background in the figure is due mainly to ip' —> 7 T°n°J/ip, 
J/ip —> pp. Fig. 1(b) shows the corresponding Monte Carlo distribution of these background 
channels, using branching fractions measured by previous experiments [7]. To remove these 
background events, we require |m pp — 3.097| > 1.5 a for 77 invariant mass (m 77 ) smaller than 
0.4 GeV/c 2 , and < 3.2—0.3m 77 for m 77 larger than 0.4 GeV/c 2 , where cr ~ 0.011 GeV/c 2 
is the pp invariant mass resolution as estimated from Monte Carlo simulation of ip' —> p'XcJ- 
After the above selection, the 77 invariant mass distributions are shown in Figs. 2(a) and 
(c), where clear 7r° and rj signals are observed. 
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FIG. 1: Scatter plots of pp invariant mass versus 77 invariant mass before removing J/%b back¬ 
ground. (a) is from data and (b) is from Monte Carlo simulated ip' —> xXcJ, Xcj 7 J/ip, J/ip PP 
{J = 0,1,2), ip' — > pJ/ip, J/ip — ► pp, and ip' —> 7r 0 7r° J/ip, J/ip —> pp events. The lines show the 
selection criterion described in the text. 


The same analysis is performed on a MC sample of 14 M inclusive ip' decays generated 
with Lundcharm [19]. It is found that the remaining backgrounds are mainly from ip' —> 
7r°7r°pp, many via resonances such as / 0 and / 2 . Some other decay channels of ip' with 
three photons such as ip' — > jXco , Xco —^ pA + ,A + —> 7r°p are also observed. Since there 
are no branching fractions available for normalization of these channels, we do not try to 
simulate all possible background channels. Instead, in our fit to the 77 mass distributions, 
we approximate the background shape by a smooth curve as predicted by the inclusive MC 
sample. 


IV. FITS TO THE 77 MASS DISTRIBUTIONS 


The 77 invariant mass distribution of candidates is not described well by a simple Gaus¬ 
sian, but by the sum of multiple Gaussians with different standard deviations, which depend 
on the momentum of the 7T° or p. The analysis of the 7T° signal is done using five different 
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FIG. 2: Invariant mass distributions of 77 from the selected ip' -► 77 pp candidate events in 
data and in Monte Carlo simulation described in Section V: (a) ppn 0 data, (b) ppir 0 Monte Carlo 
simulation, (c) ppp data, and (d) ppp Monte Carlo simulation. The curves show the best fit to the 
distributions. 

momentum bins, which are fit with multiple Gaussians for the signal and a second order 
polynomial for the background. Summing up all the fits yields the curve in Fig. 2(a), and 
the total number of ip' —> ppir 0 events is found to be 256 ± 18, with the error determined 
from the fit. 

The number of ip' —> ppp events is even more limited than ip' —> ppn°, and we do not 
do the fit in ?/ momentum bins. Instead the 77 invariant mass spectrum is fit with a single 
Gaussian for the signal plus a second-order polynomial for the background. I 11 the fit, the 
mass resolution is fixed to 14.3 MeV/c 2 , which is determined from Monte Carlo simulation 
but calibrated using the 7T° signal in the ppir 0 fit. The fit is shown in Fig. 2(c), and the 
number of events is found to be 44. 8 +g' 4 . The p mass from the fit, (552.4 ± 3.2) MeV/c 2 , 
agrees well with the world average [7]. The statistical significance of the ppp signal is 
estimated to be 6.1er by comparing the likelihoods with and without the signal in the fit. 


V. RESONANCE ANALYSIS AND EFFICIENCY 


In order to determine the selection efficiency, it is necessary to know the intermediate 
states in the decays for Monte Carlo simulation. Figs. 3(a) and (c) are the Dalitz plots for 
ip' —> ppn° and ip' —»■ ppi ] after requiring the 77 invariant mass to be consistent with a 7r° 
(0.11 GeV/c 2 < m 77 < 0.16 GeV/c 2 ) or p (0.53 GeV/c 2 < m 77 < 0.57 GeV/c 2 ). I 11 these 
figures, the requirement on the pp mass is removed to see the effect of the backgrounds 
remaining from ip' — > pJ/il >, — * 7 r° 7 r°J/' 0 , and ip' —> 7 x c j i n the lower left of the Dalitz 
plots. These two figures differ significantly from phase space. However, the data samples 
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here are too small to perform a partial wave analysis. 




FIG. 3: Dalitz plots for ip' — > ppn 0 and ip' —> pprj. (a) and (b) are for ip' — > ppir 0 data and the 
mixed Monte Carlo sample, respectively, and (c) and (d) are for ip' —> ppij data and the mixed 
Monte Carlo sample, respectively. 

In the MC simulation, iV*(1535)p + c.c. and i?(2000)7r° were used in the ppri 0 mode and 
iV*(1535)p + c.c. and i?(2000)?7 in pprj mode, where i?(2000) is a state representing the 
accumulation of the events near pp mass threshold with 

piV*(1535) : piV*(1535) : vr°i?(2000) = 2:2:1 

for the ppn 0 mode, and 

piV*(1535) : piV*(1535) : rjR( 2000) = 5:5:3 

for the pprj mode, and N*, 7r°, and rj only decay into desired final states. The MC simulations 
of the Dalitz plots are shown in Figs. 3(b) and (d). The agreement between data and MC 
simulation is reasonable. 

Fig. 4 shows the angular distributions for our selected data samples as well as the mixed 
MC samples, where 9 is the polar angle of the proton measured in the ip' rest frame, 6* and 
(p* are the polar and azimuthal angles of the anti-proton in the rest frame of pir (or prj). The 
simulations are similar to data, although not perfect. Using the same analysis as used for 
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data on the two mixed MC samples, with proper fractions of background added to the 77 
mass distributions, yields efficiencies of (14.04 ± 0.14)% for ip' —* ppir 0 and (14.00 ± 0.20)% 
for ip' —» pprj. 






FIG. 4: Angular distributions of selected ip' —> ppn 0 (top) and ip' —> ppp (bottom) events. The dots 
with error bars are data and the histograms are MC simulation (normalized to the same number 
of events). The first through third columns are cos 8, cos 9* and cp * distributions. 


VI. SYSTEMATIC ERRORS 


The systematic errors, whenever possible, are evaluated with pure data samples that are 
compared with the MC simulations. Table I lists the systematic errors from all sources. 
Adding all these errors in quadrature, the total percentage errors are 11.2% and 11.6% for 
7 T°pp and 7 ipp respectively. The detailed analyses are described in the following. 


A. Photon ID 

The fake photon multiplicity distributions and energy spectra for both data and Monte 
Carlo simulation for ip' —> pp are shown in Fig. 5. For this decay channel, the Monte 
Carlo simulates slightly less fake photons than data, while it simulates the energy spectra 
reasonably well. 

Using a toy Monte Carlo simulation, we found that for (97.10 ± 0.32)% of the cases, 
the energies of both real photons are greater than those of the fake ones from data, while 
for Monte Carlo simulation, the corresponding fraction is (97.78 ± 0.15)%. A factor of 
(0.993 ±0.004) is found between data and Monte Carlo. We do not apply a correction to the 
MC efficiency; instead, we take 1.1% as the systematic error of photon identification (ID). 






























































TABLE I: Summary of systematic errors. Numbers common to the two channels are only listed 
once. 


Source 

7 r°pp (%) 

T]pp (%) 

MC statistics 

1.0 

1.4 

Photon ID 

1.1 

Photon efficiency 


1 

7r°(r/) reconstruction 

2.0 

Tracking and particle ID 

2.6 

2.8 

Fit to mass spectrum 

4.5 

5.3 

Decay dynamics 

5.9 

Kinematic fit 



Number of ip' 


1 

Trigger efficiency 

0.5 

“*• 77] 

0.0 

0.7 

Total Systematic error 

11.2 

11.6 


600 

500 

$400 

§300 
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FIG. 5: Fake photon multiplicity distributions (left) and fake photon energy spectra (right) for 
ip' —> pp data (dots with error bars) and Monte Carlo simulation (histogram). 


B. Photon detection efficiency 


The simulation of the photon detection efficiency is studied using J/ip —> 7r + 7r _ 7r° events 
with one photon missing in the kinematic fit and examining the detector response in the 
missing photon direction [20]. The Monte Carlo simulates the detection efficiency of data 
within 2% for each photon in the full energy range. Since we have two photons, 4% is taken 
as the systematic error of the photon detection efficiency. 


9 



















C. 7T° and p reconstruction 


The 7 r° reconstruction efficiency is studied by comparing the opening angle between the 
two photons between data and MC simulation in different 7T° momentum ranges using J/ip —> 
ppir 0 and J / pi —> 7r + 7r _ 7r° samples selected from BES J/'ip data. Fig. 6 shows the comparison 
for 7r° momentum between 0.5 GeV/c and 0.6 GeV/c, the agreement at small opening 
angle shows the simulation of the 7T° reconstruction efficiency is good. By reweighting the 
difference between data and MC simulation in all momentum bins with the 7T° momentum 
spectrum in xp' —> ppir 0 , the overall correction factor to the MC simulation is determined to 
be (98.8±0.8)%, and 2.0% is then taken as the systematic error due to the 7T° reconstruction. 



FIG. 6 : Comparison of the 77 opening angle (a) distributions for 7 r° from J/ip with 

momentum between 0.5 GeV/c and 0.6 GeV/c. Dots with error bars are data, and the histogram 
is Monte Carlo simulation. The distributions are normalized to the number of events with a > 28°. 

The angle between the two photons emitted from the p is generally much greater than 
that from the 7 T°. As a conservative estimation, the uncertainty for p reconstruction is also 
taken to be 2 . 0 %. 


D. MDC tracking and particle ID efficiency 


The efficiencies for protons and antiprotons that enter the detector being reconstructed 
and identified are measured using samples of J/ip — > tt + tt ~pp and ip' — » tt + n~ J / ip, J /ip — > pp 
events, which are selected using kinematic fit and particle ID for three tracks, allowing one 
proton or antiproton at a time to be missing in the fit. The efficiency is determined by how 
often a proton (or antiproton) is found in the direction of the missing track, it varies from 
80% to around 95% with increasing momentum, and the MC simulates data rather well, 
except for proton or antiproton momenta less than 0.5 GeV/c, where the nuclear interaction 
cross section of particles with the detector material is very large. 

The net difference between data and MC simulation is found to be 1.001 ± 0.025 for 
ip' —> ppir 0 , and 1.010 ± 0.018 for ip' —> ppp. The errors together with the differences from 
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unity will be considered as systematic errors, that is, 2 . 6 % and 2 . 8 % for i\)' —> ppn° and 
0 ' —>• pprj, respectively. 


E. Fit range and background shape 


The background shape in fitting the 77 mass distributions is changed from a second order 
polynomial to a first order one, and the fit range is changed, to determine the uncertainties 
due to the fitting for the ppn° and pprj channels. Different ways for choosing the ir° momen¬ 
tum bins or fitting the 7r° signal without binning yields differences in the branching fraction 
less than 3% for the ppn° channel. Adding all these in quadrature, 4.5% and 5.3% are taken 
as the systematic error due to the fit. 


F. Decay dynamics 


Table If shows efficiencies determined with different MC samples; different decay dy¬ 
namics result in different efficiencies. While the mixed Monte Carlo samples with N* and 
possible pp intermediate states are used in the determination of final selection efficiencies, 
the differences between the mixed samples and the phase space generator are taken as sys¬ 
tematic errors due to the lack of the precise knowledge of the decay dynamics to be used 
in the MC generator. The differences are found to be 2.1% for 0 ' —> ppn° and 5.9% for 
0 ' —>• pprj. The larger difference (5.9%) will be taken as the systematic error for both due to 
the uncertainty of the generator for both channels. It should be noted that the differences 
between these two MC samples in the angular distributions are large compared with those 
observed in Fig. 4; thus the errors quoted cover the differences in the angular distribution 
simulation also. 


TABLE II: Efficiencies determined with different MC samples. 


Channel 

ppn 0 (%) 

pprj (%) 

Only A* (1535) 

13.04 

12.71 

Only f?(2000) 

18.43 

18.40 

Mixed sample 

14.04 ±0.14 

14.00 ±0.20 

Phase space 

14.37 

14.83 


G. Other systematic errors 

The uncertainty due to the kinematic fit is extensively studied using many channels 
which can be selected cleanly without using a kinematic fit [21-24], It is found that the MC 
simulates the kinematic fit efficiency at the 5% level for almost all the channels tested. We 
take 5% as the systematic error due to the kinematic fit. 
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The results reported here are based on a data sample corresponding to a total number of 
ip' decays, TV,//, of (14.0 ±0.6) x 10 6 , as determined from inclusive hadronic events [15]. The 
uncertainty of the number of ip' events, 4%, is determined from the uncertainty in selecting 
the inclusive hadrons. 

The trigger efficiency is around 100% with an uncertainty of 0.5%, as estimated from 
Bhabha and e + e _ —■> /i + /V events. The systematic errors on the branching fractions used 
are obtained from the Particle Data Group (PDG) [7] tables directly. 


VII. CONCLUSION AND DISCUSSION 


The branching fractions of ip' —> ppn° and ip' —»• pprj are calculated using 

n$/e 


B(ip' —> ppn°) = 
B(ip' —> pprj) = 


■ B( 7 T° -> 77) ’ 


n° b s /e 


Ny ■ B{j] —> 77 )' 


Using numbers listed in Table. Ill, we obtain 

B(ip' —> pp7r°) = (13.2 ± 1.0 ± 1.5) x 10" 5 , 

B(ip' -> pprj) = (5.8 ± 1.1 ± 0.7) x 10~ 5 , 

where the first errors are statistical and the second are systematic. The measured ip' —► pp7r° 
branching fraction agrees with Mark-II within errors [14]. 


TABLE III: Numbers used in the branching fraction calculation and final results. 


quantity 

ppn 0 

pprj 


n °bs 

256 ± 18 

44.8+U 


e{%) 

14.04 ±0.14 

14.00 ± 0.20 


N^( 10 6 ) 

14.0 ±0.6 


°(h)-77)(%) 

98.80 ± 0.03 

39.43 ± 0.26 

B(ip' - 

-* ppir°(ri)) (10 -5 ) 

13.2 ± 1.0 ± 1.5 

5.8 ± 1.1 ±0.7 


Comparing the branching fractions with those of J/ip decays from the PDG [7], we find 
that for ip' decays, B(ppn°) > B(pprj), while for J/ip decays, B(ppn (> ) < B(ppij), and 


Qppn 0 

Qppn 


B(ip' — > ppn°) 
B(J/ip —» pp-ir 0 ' 
B(ip' —> pprj) 
B(J/ip —» pprj) 


= (13-2 ± 1-0 ± 1-5) x 10- = (121 ± 19)% 
(1.09 ± 0.09) x lO " 3 V ; ’ 

(5.8 ± 1.1 ± 0.7) x 10~ 5 / 9 q 1 n 7 w 


While Q P p n 0 agrees well with the so-called 12% rule, Qppp seems to be suppressed significantly. 
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Fig. 7 shows the pp invariant mass distributions of the selected ppn° and ppp events shown 
in Fig. 3, together with the expected background estimated from 7r° or p mass sidebands 
(0.075-0.100 and 0.170-0.195 GeV/c 2 for vr° and 0.49-0.51 and 0.59-0.61 GeV/c 2 for p). 
There are indications of some enhancement around 2 GeV/c 2 in both channels. Fitting the 
enhancement with an S-wave Breit-Wigner and a linear background, with a mass dependent 
efficiency correction, yields a mass around 2.00 GeV/c 2 in the ppn 0 mode and 2.06 GeV/c 2 
in ppp , with the width in both channels around 30-80 MeV/c 2 , and significance around 2.7cr. 
Fitting with a P-wave Breit-Wigner results in slightly lower masses and similar significance. 
The nature of the enhancements is not clear, and the statistics are too low to allow a detailed 
study. The enhancements in the two channels cannot be the same since they have different 
isospin. 




FIG. 7: pp invariant mass distributions of selected (a) ppir 0 and (b) ppp events. The blank his¬ 
tograms are selected signal events, and the shaded histograms are events from n r° or p mass side¬ 
bands. The dashed histograms are predictions of phase space with S-wave pp (not normalized). 

Fig. 8 shows projections of Dalitz plots in pir (or pp) invariant mass after removing 
ifj' —> J/^ + X backgrounds and the possible pp mass threshold enhancements. There is a 
faint accumulation of events in the pit invariant mass spectrum at around 2065 MeV/c 2 , but 
it is not statistically significant. The enhancement between 1.4 and 1.7 GeV/c 2 may come 
from 7V*(1440), 1V*(1520), 1V*(1535), etc. We do not attempt a partial wave analysis due to 
the limited statistics. There is a clear enhancement with pp mass at (1549 ± 13) MeV/c 2 , 
which is possibly the IV* (1535). 


VIII. SUMMARY 


pp7r° and ppp signals are observed in / decays, and the corresponding branching fractions 
are determined. For ijj' —> ppn°, the errors are much smaller than those of the previous 
measurement by Mark-II [14], and for ifr' — > ppp , it is the first observation. There is no 
clear IV* (2065) peak in the ppir 0 mode, but there is some weak evidence for pp threshold 
enhancements in both channels. 
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pn° mass (GeV/c 2 ) pn 0 mass (GeV/c 2 ) p7t°/p7i° mass (GeV/c 2 ) 
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pr) mass (GeV/c 2 ) p-q mass (GeV/c 2 ) pr|/pri mass (GeV/c 2 ) 


FIG. 8: Projections of Dalitz plots in jm{prj) invariant mass after removing ip' —> r]J !> and the 
possible i?(2000). (a) and (b) are m pn o and m pn o in %/ —> ppir 0 ; (c) is the sum of (a) and (b); (d) 
and (e) are m pv and m PV in ip' —> ppp-, (f) is the sum of (d) and (e). 

Acknowledgments 

The BES collaboration thanks the staff of BEPC for their hard efforts. This work is 
supported in part by the National Natural Science Foundation of China under contracts 
Nos. 10491300, 10225524, 10225525, the Chinese Academy of Sciences under contract No. 
KJ 95T-03, the 100 Talents Program of CAS under Contract Nos. U-ll, U-24, U-25, and 
the Knowledge Innovation Project of CAS under Contract Nos. U-602, U-34 (IHEP); by the 
National Natural Science Foundation of China under Contract No. 10175060 (USTC), and 
No. 10225522 (Tsinghua University); and by the Department of Energy under Contract No. 
DE-FG02-04ER41291 (University of Hawaii). 


[1] J. L. Rosner, Phys. Rev. D 64, 094002 (2001). 

[2] P. Wang, C. Z. Yuan and X. H. Mo, Phys. Rev. D 70, 114014 (2004). 

[3] N. Isgur, G. Karl, Phys. Rev. D 18, 4187 (1978). 

[4] N. Isgur, G. Karl, Phys. Rev. D 19 , 2653 (1979). 

[5] K. F. Liu, C. W.Wong, Phys. Rev. D 28 , 170 (1983). 

[6] S. Capstick, N. Isgur, Phys. Rev. D 34, 2809 (1983). 

[7] Particle Data Group, S. Eidelman et al., Phys. Lett. B 592 , 1 (2004). 

[8] B. S. Zou, Nucl. Phys. A 675 , 167C (2000). 

[9] BES Collaboration, J. Z. Bai et al., Phys. Lett. B 510 , 75 (2001). 

[10] BES Collaboration, M. Ablikim et al., hep-ex/0405030. 


14 








[11] S. Capstick and W. Roberts, Phys. Rev. D 47, 1994 (1993). 

[12] S. Capstick and W. Roberts, Prog. Part. Nucl. Phys. 45, S241 (2000). 

[13] C. H. Chang and H. R. Pang, hep-ph/0407188. 

[14] Mark-II Collaboration, M. E. B. Franklin, et ai, Phys. Rev. Lett. 51 , 963 (1983). 

[15] X. H. Mo et al, HEP&NP 28, 455 (2004). 

[16] BES Collaboration, J. Z. Bai et al., Nucl. Instr. Meth. A 344, 319 (1994). 

[17] BES Collaboration, J. Z. Bai et al., Nucl. Instr. Meth. A 458, 627 (2001). 

[18] BES Collaboration, M. Ablikim et al., physics/0503001. 

[19] J. C. Chen et ai, Phys. Rev. D 62, 034003 (2000). 

[20] S. M. Li et al, HEP&NP 28, 859 (2004). 

[21] BES Collaboration, J. Z. Bai et al., Phys. Rev. D 69 , 092001 (2004). 

[22] BES Collaboration, M. Ablikim et al., hep-ex/0408047. 

[23] BES Collaboration, J. Z. Bai et al., Phys. Rev. D 70, 012005 (2004). 

[24] BES Collaboration, J. Z. Bai et al., Phys. Rev. D 69 , 072001 (2004). 


15 



